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Abstract 

We present a new solution of the asymmetric two-matrix model in the large N limit 
which only involves a saddle point analysis. The model can be interpreted as Ising in 
the presence of a magnetic field, on random dynamical lattices with the topology of the 
sphere (resp. the disk) for closed (resp. open) surfaces; we elaborate on the resulting phase 
diagram. The method can be equally well applied to a more general (Q + l)-matrix model 
which represents the dilute Potts model on random dynamical lattices. We discuss in 
particular duality of boundary conditions for open random surfaces. 
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1. Introduction 



The study of various multi-matrix models in the large iV limit is motivated by their 
interpretation as statistical lattice models on random surfaces. The one-matrix model [|1| 
already describes the summation over random surfaces, but to put "matter" on the surface, 
several matrices are requiredB. The simplest such model is the two-matrix model, which 
has the following partition function ||,[5| 



iVtr 



\{A 2 + B 2 ) + ^fA 3 + + ±AB 



Z(a ,A>,7) = jj dAdBe 1 7 J (1-1) 

where A and B are N x N hermitean matrices. In the large iV limit, this generates 
triangulated surfaces with the spherical topology, on which spins live (the two matrices A 
and B correspond to spins up and down) [[7| , thus reproducing the Ising model on random 
surfaces. Schematically, ao and /3q play the roles of both magnetic field H and cosmological 
constant (ao/ flo = exp(2H/T)), while 7 is related to the temperature T by 7 = exp(l/T). 

This model has been solved using orthogonal polynomials |J, but, strangely, the 
simplest tool available, the saddle point method, which works so well for the one-matrix 
model, has not been used. It is commonly assumed that this method does not work in 
the two-matrix model, the usual argument being the following: after diagonalization of A 
and B and use of the Itzykson-Zuber-Harish Chandra formula we are left with an 
integral over the eigenvalues: 



Z(a ,A),7) = J J n da * d6 * A N A [^] det 



e 7 3 



c 



(1.2) 

where A(-) is the Van der Monde determinant. Then, one uses symmetry of permutation 
of the eigenvalues to reduce this expression to: 

7( R \ /YtTh Ah Al lAfM N T,i\-h(« 2 i + + 

Z(o ,Po,7)= // [[daidbiA[ai\A[bi\e L 7 

(1.3) 

At this stage, a large N saddle point analysis shows that there is a continuous infinity of 
saddle points and it is very difficult to derive anything from it. Of course, the problem 
stems from the transformation of ( |1.2|) to (|0|): since we have broken the symmetry of 



1 Except the Kazakov multicritical points || which correspond to non-unitary matter. 
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permutation of the eigenvalues, each ordering a CT (i) < • • • < o,ct{N) of the eigenvalues leads 
to a different saddle point, so that we have AH saddle points, which causes trouble as 
A" — > oo. The problem does not exist at the level of ( |1.2| ) and in fact, as we shall show, 
there is a well-defined saddle point to it; we shall then explain how to determine it, which 
will result in a very compact and elegant way of expressing the resolvent (s) of the model. 
A similar result could be obtained for the two-matrix model with quartic vertices, but we 
shall not choose to do so here. 

We shall then generalize our results to the dilute Q-states Potts model on random 
surfaces, which is defined by (see [E21 for a somewhat similar definition on a flat lattice) 



Q JVtr 



Zq(«o,A),7)= f[dAf[dB q e 

J J 0=1 



(1.4) 

This model describes the following statistical model on random surfaces: each vertex of the 
surface is either unoccupied (represented by matrix A) or occupied by a spin in Q possible 
states (matrices B q ), with the rule that adjacent vertices cannot be occupied by spins in 
different states. Again, «q and (3q are cosmological constants and control the dilution (i.e. 
density of unoccupied sites) , and 7 is related to the inverse temperature. The two-matrix 
model ( |1.1|) is the particular case Q = 1. 

From now on we shall redefine the couplings of the model and rescale the fields so 
that the partition function can be rewritten 



Z^ft^/cUe^I-^ + f^] (J dBe N «l-l B2 + % B3 + AB \\ Q (1 .5) 
where a = ao7 

2, p = /3Q7 3 / 2 . The main physical quantities of the model are the 
resolvents, which are defined by 

/ \ / 1 1 
u A {a) = ( — tr 



. N a-A. 
^ (6) = (iV tr ^TB 

where the large A^ limit is implied and a, b are complex numbers. u^(a) and uJB(b) 
are generating functions of averages of the form (tr A n ) and (trB n ), but they are also 
important from the diagrammatic point of view: they correspond to sums over connected 
surfaces with a boundary ("loop functions"), the parameters a or b playing the role of 
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boundary cosmo logical constant. In the large N limit, these surfaces have the topology 
of a disk. The difference between ua{o) and us(b) lies in the boundary conditions: for 
u)a{o) (resp. u>B(b)), there are only matrices A (resp. B) on the boundary. The large n 
asymptotics of tr A n and tr B n (i.e. surfaces with large boundary) are dominated by the 
singularities of the corresponding resolvent. These singularities are also relevant for the 
physics in the bulk: if g is the exponent of this singularity, then the central charge c of the 
critical model is given by [|18| , |T9 l 



c = l-6{ y /g-l/^g) 2 (1-7) 

These resolvents will therefore play a central role in our analysis, and our goal will be to 
find exact expressions for them. Let us remark that when the dilution is turned off (a = 0), 
one can perform the gaussian integration over A and we are brought back to the usual 



Q-states Potts model on random surfaces [|15|. cus(6) is then the standard resolvent: it 



corresponds to boundary conditions of the Potts model where all the spins on the boundary 
are in a given state. ua{o) is not a natural resolvent, because the sites on the boundary are 
unoccupied, whereas they cannot be so in the bulk. However, once dilution it turned on, the 
two resolvents ua(o) and o;b(6) should be treated on equal footing; they will correspond, 
as we shall explain later, to boundary conditions which are "dual" to each other. 

The plan of the article therefore goes as follows: first we shall analyze in section 2 
the integral over B which is common to the models Ql.5| ) for arbitrary Q, then discuss the 
cases Q = 1, 2 and 3, 4 (though in ( |1.5| ) the parameter Q can take arbitrary real values, 
for simplicity we only consider here integer values) in sections 3, 4, 5 and finally conclude 
in section 6. 



2. The external field problem 



Let us first consider part of the model only: if the matrix A is held fixed in (1.5), we 



are left with the problem of one matrix in an external field [|13l , |14 

E(A)=fdBe Nt ^- V (B) + AB] 



(2.1) 



where V is a polynomial potential. We shall show how the saddle point equations expressed 
JTT|] allow indeed to calculate this matrix integral, and how in the case of the cubic 



m 



potential V(b) = — § b 3 + ^6 2 , they in fact reproduce the solution |14[] without the use of 
any partial differential equations. 
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By [/(iV)-invariance only depends on the eigenvalues of A: 

detfe^i] 



II ^[fc] " A[oj e-^ ^ W (2.2) 



We shall now show how to calculate in the large N limit the logarithmic derivatives of S 
with respect to a, (S then follows by simple integration). 

We assume that the density of eigenvalues of A becomes smooth in the large iV limit, 
and denote it PAifl); for the sake of simplicity only, its support will be taken to be of the 
form of a single interval [ai, 02]. It is related to the resolvent of A by 

UA{a) = / — (2.3) 

which is an analytic function everywhere except on the support of A, where it has a cut 
(leading, if pa{&) is smooth, to other sheets) which we call the physical cut. 



We can now write 11 



1 d 

T?^— !og s k] = b(a,i) - u A {a,i) (2.4) 

JS OCLi 

where 6(a) is an analytic (multi- valued) function which has the same physical cut as ua(o)- 
The particular sheet corresponding to the value of 6(a) in ( |2.4j) is called the physical sheet. 
We see that 6(a) is the quantity we need to compute. 

In order to do so, we must write down the saddle point equation of ( |2.1|) , which we 
shall do carefully. First we introduce symmetrically the resolvent wb(6) of B, with a 
physical cut [61, 62] (i.e. the eigenvalues of B fill the interval between 61 and 62): 

^ &b'p B {b') 



uj B {b) 



fcx b-b> 



Taking the logarithmic derivative with respect to the bi results in the appearance of the 
function a(6), which has the same cut as the resolvent of B, and is the functional 

inverse of 6(a) (see appendix 1 of § jj). The saddle point equations then read: 



\{u B {b + iO) + uj B (b - iO)) + ha(b + iO) + a(b - iO)) = ^(6) 6 e [b u b 2 ] 

Using the fact that a(6) and u) B (b) have the same cut, this equation can be analytically 
continued: 

u B (b) + a*(b) = V'(b) (2.5) 
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where a* (6) is the value of a(b) on the other side of the physical cut of B. Before going 
further, let us give the physical significance of a(b). Formally, the saddle point equation 



(|2.5|) can be written down: ^ <5 P tf(b) = ^' wnere S is the action, and the derivative d/d& 



takes care of the normalization condition J dbpsib) = 1. Now SS/dpsib) is nothing but 
the effective potential V e s(b) for the eigenvalues bi, which must be the sum of three terms: 
the potential V(b), the effective interaction among the 6j, and the effective interaction 
between the b{ and the Oj. Indeed, we find explicitly that 

{const b £ [&i, 62] 

,6 
y(6)- / [w B (6) + a(6)]d6 b<£[b x M 

(note that it is a(b) which appears and not a* (6) because we are outside the cut). Therefore 
a(b) can be interpreted as the derivative of the effective potential for the action of the at 
on the &Ji. 

Let us show now how to handle the equation (|2.5|) : it can be rewritten 



u B (b(a))+a = V'(b(a)) (2.6) 

where it is understood that ( |2.6| ) is a relation between multi-valued functions. To use this 
equation, we now assume that V is a polynomial of degree d, and look at the behavior of 
b(a) as a — > 00. Under minimal assumptions on the analytic structure of b(a) (see below), 
on all cuts connected to the physical sheet without crossing the physical cut, b(a) — > 00 
as a — > 00. Therefore b(a) d ~ 1 ~ a as a — > 00, which leads us to the "minimal" conjecture 
for the analytic structure of b(a) which is depicted by figure fig. 1 (a similar structure was 
found for large iV characters in |T7[]). All the cuts except the physical cut go to infinity, 
and we assume that they do not cross the physical cuti. 



2 Symmetrically, b(a) will be the derivative of the effective potential for the action of the bi on 
the a,i. 



This corresponds to the strong coupling phase of [13|. In the present context, it simply means 
we remain below the continuum limit surface. 



physical cut 



d-1 sheets 



^ ^—^1 ! " 

















physical sheet 



Fig. 1: Analytic structure of b(a) in the external field problem with a polynomial 
potential of degree d. 

At this point we have complete knowledge of the analytic behavior of 6(a) (the physical 
cut is determined by the fact that it is identical to that of uja(cl), which is known by 
definition: b(a — iO) — b(a + iO) — 27iipA(a)), and this is enough to compute it. 

More explicitly, by a change of variables of the form: a = P(z), where P is a polyno- 
mial of degree d — 1 whose critical values are the branching points, we can remove all the 
cuts at infinity; then b(z) has a single cut, the physical cut, and it can be expressed as: 



b(z) = c\z + c 2 + 



z ^ dz'p A (a(z')) 

(ai) 



Z — Z 



(2.7) 



where the constants c\ and c 2 are easily determined by asymptotics at infinity. 

The simplest non-trivial case is the cubic case, in which we find only two sheets (the 
physical sheet and an extra one) if we do not cross the physical cut, and the equation ( |2~7| ) 
reduces to the known solution WM. Since this is the case that is of interest to us, let us 
carry out explicitly the procedure outlined above. From Eq. ( |2.6| ) with V(b) = — f 6 3 + ^6 2 , 
we obtain that there are two sheets b±(a) (the physical sheet being by definition 6+ (a)) 
on which as a — > oo, 



-1/2 



(2.! 



where we have used w_b(6) = 1/6 + 0(l/6 2 ).@ After the change of variables z 2 = 
(where a3 is the branch point of the semi-infinite cut), we find that 



4a 3 (3 1 

z 



?3/2 



2z 2 



+ 0(l/z 3 ) 



Had we used an expansion o;_b(6) = ^^-0 B n /b n+1 where B n = (jjtrB n ^, we would have 
obtained an expansion at arbitrary order of b±(a) in terms of the B n . 
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In particular the constants in (|2.7|) are given by c\ = ft -1 / 2 and C2 = -^s- 

Note however that in the next paragraphs, we shall not need to make the explicit 
change of variables a — > z. It is known that this maps the saddle point equations of the 



Q-states Potts model onto those of the 0(n) model |16[], but this correspondence has no 
direct physical meaning (the phase diagrams are different, the critical models do not have 
the same central charge, etc) and is more misleading than useful for our purposes. Instead, 
we shall directly consider 6(a) in its normal parametrization and use the information on 
its analytic structure discussed above. 



3. The Q = 1 dilute Potts model (Ising with magnetic field) 

Let us now apply the results of the previous paragraph to the two-matrix model. The 
philosophy will be the same in this case as in the more general Q-states Potts model with 
< Q < 4: we have two saddle point equations, one coming from A, the other one from 
B. These equations involve functions a(6) and b(a) which satisfy a functional inversion 
relation; therefore, by inverting one of them, we obtain two relations for 6(a). Recombining 
them, we deduce a polynomial equation satisfied by 6(a). 

The saddle point equation of the eigenvalues of B has already been analyzed in the 
previous section; as we have shown, it gives the behavior of 6(a) at infinity and fixes its 
analytic structure (Eq. (|2.8|)). We have found two sheets, the physical sheet 6_|_(a), and 
6_(a) which is connected to it by a semi-infinite cut. 

Remembering that the full partition function takes the form 

Zl («,/3, 7 ) = // n^dkAMAfe] det [e^e^ N (a? + ^ + ^ + 

in terms of the eigenvalues of A and B, we now write the analytically continued saddle 
point equation for the eigenvalues of A: 

u>a(ci) + 6* (a) — 7a + aa 2 = (3-1) 

6* (a) is a third sheet of 6(a), connected to the physical sheet 6 + (a) by the physical cut, 
and which according to (|3.1| ), has no other cut than the physical cut. This means that 6(a) 
has exactly three sheets, and we therefore make the Ansatz that 6(a) is the solution of a 
third degree equation in a. The coefficients of this polynomial are symmetric functions of 
the different sheets; for example, using the elementary identity 6+ (a) +6_(a) = uja(cl) + -k 
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(coming from the analytic structure of b(a)) and (|3.1|) , we find that b + (a) +6_(a) +6* (a) = 
— aa 2 + 7 a + 5. For other symmetric functions we must use the finite expansion (|2.8|) , 
which implies that not all of the coefficients of the polynomial can be found explicitly. We 
find that there are three remaining unknowns called x, y, z (which can be reexpressed in 
terms of the three first moments of 6, cf footnote 4 ). The final equation can be cast in the 
nicely symmetric form: 

a(3 a 2 b 2 + a a 3 -<yaa 2 b-<y[3b 2 a+/3 b 3 -7 a 2 + (<y 2 + l-a(3) ab-<y b 2 + x a+yb+z = (3.2) 

The constants x, y, z are determined by imposing that a(b) and b(a) have the appropriate 
analytic structure: when solving ( |3.2| ) for 5, the discriminant is a 9 th degree polynomial in 
a, and we must impose that 3 zeroes are double zeroes in order to have only three branch 
points!. This gives three algebraic equations for x, y, and z which are therefore functions 
of a, f3, 7. 

One can show that, in the absence of magnetic field (a = 0), the equation (|3.2j ) 
is equivalent to the equation for the resolvent found in |8|,|9"l,|lCi| using more complicated 
methods; however, (|3.2| ) displays explicitly the Z2 symmetry a <-> b, a «-> f3, whereas it is 
not obvious in P,BL[TO]. This explicit symmetry is due to the fact that we are not trying 
to write down an equation for ua(ci) directly (or ws(6)), but rather for 6(a), which differs 
by a polynomial part —7a + aa 2 . 

Figure fig. 2 shows the resulting topology for the phase diagram of the model. We 
shall now explain it by briefly analyzing Eq. (|3.2|) in two particular cases. 




1/7 



Po 



Fig. 2: Schematic phase diagram of the Ising model with magnetic field on random 
surfaces. PT is the zero magnetic field low temperature phase transition line; T is 
the critical point of the model; CT (resp. DT) is a boundary phase transition line 
for b(a) (resp. a(b)). 



5 Note that this constraint is much stronger than the one given in |] 
which in fact would not be enough to fix the three unknowns. 



' "one cut hypothesis" ' 
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If we first remove the dilution (a = 0, which in the Ising language, corresponds to 
an infinite magnetic field), this model is simply the one-matrix model, and we expect no 
physics at all; in fact, integrating over A shows that there is only one parameter in the 
model, /?/(7 — I/7) 3 / 2 . The continuum limit, which is attained for 

P 2 1 



(7-1/7) 3 12^3 

corresponds to a c = theory, and there seems to be no critical point. However, this is 
wrong, because even though the bulk theory is pure gravity with a coupling 0/(^y — 1/7) 3 / 2 , 
6(a) represents a non-trivial loop function, and the corresponding boundary operator de- 
pends explicitly on 7. Indeed we find that the standard resolvent wb(6), or equivalently 
a (6) (cf Eq. (|2.5|)) always has the singularity 



a — a reg ~ 



(6-M 3/2 



of the usual pure gravity loop function, but that 6(a) (or ua(o)) undergoes a phase transi- 
tion! This can be seen in the behavior of the branch points of 6(a) as one varies 7 (figure 
fig. 3). There is a critical point (point C of figure 2) 



7c 2 = 1 + 2V3 ai = 2 7c 

where the two cuts merge, and the result is that 



3 





((a 


-a 2 , 


(3/2 


7 < 7c 


6 - 6 reg ~ | 




-a 2 ] 


,2/3 


7 = 7c 




[(a 


-a 2 ] 


jl/2 


7 > 7c 



(3.3) 



For 7 > 7c there is again a (a — a^) 3 ^ 2 singularity, but as does not belong to the physical 
cut. 

113 14 13 1 

• • • • • • • • 

a ! a 2 a 3 a i a 2= a 3 a i S S 

Y>y c 'fy, y<y c 

Fig. 3: Cut structure of b{a) for a = (infinite magnetic field) on the continuum 
limit surface. The numbers represent the multiplicity of the zeroes of the discrimi- 
nant. Note that the two cuts never intersect each other (the choice of cuts is dictated 
by analytic continuation from the gaussian model). 
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Let us insist that the corresponding theory on the sphere (i.e. the free energy of 
the matrix model) only depends on the combination /3/(7 ~ I/7) 3 / 2 and is always pure 
gravity; but there is nonetheless a phase transition, at the boundary of the surface, which 
is signalled by a change of analytic behavior of b(a) from 7 < 7 C (high temperature) to 
7 > 7c (low temperature) - or equivalently, a change in the asymptotics of (tr A n ) as 
n — > 00. The situation might seem trivial at the point C where the integration over A 
is gaussian, but in fact the phase transition occurs on a whole "critical line" (line CT of 
figure 2), and can be interpreted as follows: inside the region PCT, where the temperature 
is low and the magnetic field favors spins B, the boundary of the random surface is made 
of spins A, whereas the bulk of the surface is mostly made of spins B, so that the boundary 
tries to avoid touching the rest of the surface by "collapsing on itself" (which results in a 
change of its Hausdorff dimension). This is a quantum gravity phenomenon which has no 
precise analogue on a flat lattice. Symmetrically, a (6) undergoes the same phase transition 
on the line DT; in fact, it is clear that for any Q < 4, the point D exists (its position being, 
up to a trivial rescaling of 7, independent of Q) and is the endpoint of a boundary phase 
transition line, with the same critical properties as for Q = 1; so that we shall not mention 
this line in the subsequent discussion of Q = 2, 3, 4. 

Second, let us discuss the physics on the zero magnetic field line a = (3, on which 
the Ising critical point (point T of figure 2) is. Note that the symmetry of the equation 
imposes then that x = y, so that there are only two unknown constants left in (|3.2|) (as 
opposed to the three found in [|10| ) . A very nice picture emerges out of the patterns of the 
zeros of the 9 th de gree discriminant (figure fig. 4), which in particular illustrates the Z2 
spontaneous symmetry breaking in the low temperature phase. 



1 3 2 3 1 8 1 3 1 

• • • • • • • • • 



6 6 

y>y T ^ y<y t 



• 2 



Fig. 4: Cut structure of b(a) for a = (3 (zero magnetic field) on the continuum 
limit surface. The numbers represent the multiplicity of the zeroes of the discrim- 
inant. The exchange of A and B, represented by the arrows, acts non-trivially 
on the branch points in the low temperature phase (7 > 7 T ), but not in the high 
temperature phase (7 < 7 T ). 
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At the critical point 

a T =(3 T = VlO = 3.16228 ... 7 T = 2V7 + 1 = 6.29150 . . . 

the resolvent develops a singularity (a* = = as) 

b-b reg ~(a-a^ 3 (3.4) 

which is characteristic of a c = 1/2 theory coupled to gravity. Let us note that since the 
resolvent must have a cubic singularity, three is the minimum number of sheets required 
of 6(a); therefore the cubic two-matrix model is the simplest possible realization of the 
c = 1/2 theory coupled to gravity. 

Since the Ising model on random surfaces has already been studied in great detail, we 
shall not elaborate any further and go back to the general case with arbitrary Q. Noting 
that the partition function ( |1.5| ) takes the form 

Z a( a, A7 ) = /d^e Artr [-^ 2 + f-4 3 ] S ( J 4)« 

where S(A) is defined by (|2.1[) with the usual cubic potential, and using (|2.4| ), we can 
immediately write down the saddle point equations for the eigenvalues of A; after analytic 
continuation, we obtain: 

(2 - Q)u A {a) + K(a) + (Q - l)b(a) =ja- aa 2 (3.5) 

where 6(a) = 6+ (a) is evaluated on the physical sheet. For Q = 2 and 3 the resolution is 
very similar to the Q = 1 case, and the resulting phase diagrams are of the same type as 
2, except that the line CT becomes a real critical line (for the bulk theory); we shall now 
give the main results, skipping the technical details. 



4. Q = 2 and 3 dilute Potts models 

The matrix model corresponding to the Q = 2 dilute Potts model on random surfaces 
is nothing but the Z2 symmetric three matrix chain ||, A playing the role of the central 
matrix and B of the two matrices at the ends of the chain. If we remove the dilution, 
i.e. set a = 0, we are back to the Ising model without magnetic field. However, adding 



dilution allows us to reach the tricritical point of the Ising model [EH 
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The saddle point equation ( |3.5|) reads for Q = 2: 



6*(a) + 6(a) = 7a - aa 2 (4.1) 

which shows that 6* (a) and b(a) must have the same cuts; since we know from section 2 
that 6(a) = 6+(o) is connected by a semi-infinite cut to 6- (a), 6* (a) = 6*+ (a) must likewise 
be connected to a fourth sheet 6*_(a). We therefore assume that 6(a) satisfies a fourth 
degree equation and compute symmetric functions of its 4 sheets. Using once again the 
expansion ( |2.8|) , the saddle point equation ( [4.1] ) and various other relations coming from 
the analytic structure of 6(a), we find the following algebraic equation: 

(3 6 4 - 2/3 ( - a a 2 + 7 a) 6 3 + (a 2 (3 a 4 - 2a(3-f a 3 + 7(^7 - a) a 2 + {-a(3 + 2 + 7 2 ) a + ■ ■ • ) 6 2 
- (7a 2 a 4 + a(a/3 - 2(1 + 7 2 )) a 3 + • • • ) 6 + (a 2 a 5 - 2a 7 a 4 + • • • ) = 

(4.2) 

where • • • means that there are some lower order terms whose coefficients must be fixed by 
the analytic structure. Note that as an equation for a, ( |4.2|) is a fifth degree equation. 

Fixing the unknown coefficients in ( fl.2| ) allows us to easily find the critical line, which 
is characterized by the collision of the physical cut with the semi-infinite cut. Its two 
endpoints are the zero dilution (a = 0) critical point: 

7c 2 = 2(1 + V7) (3 2 c = IO7- 3 

where the loop functions display the same behavior that is found along the whole critical 
line, namely 

6-6*~(a-a0 3/4 

(4-3) 

a - a, ~ (6 - 6,) 4/3 
and are characteristic of a c = 1/2 theory; and the tricritical point:i 

a T = 2.83045 .. . (3 r = 3.09138 .. . 7 T = 6.23472 .. . 

where the corresponding singularities of the loop functions are 

6 - 6 reg ~ (a - a*) 5/4 

(4.4) 

a - a reg ~ (6 - 6*) 5/4 

6 Only the numerical values are given, the exact values being fairly cumbersome. This remark 
also applies to the cases Q = 3 and 4. 
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As expected this corresponds to the central charge c = 7/10 of tricritical Ising. 

The 3-states dilute Potts model is the first in which the corresponding matrix model 
(5 = 3 does not have the form of a linear chain; in particular, it is not solvable via 
orthogonal polynomials. We start once more from the saddle point equation ( fe.5| ). The 
discussion of the analytic structure of 6(a) in this case becomes a little more involved, 
and let us simply state the conclusion of this analysis, which is that 6(a) must have 6 
sheets, as shown on figure fig. 5. The reader can check that this analytic structure, and 
in particular the discontinuities shown, are compatible with the saddle point equations. 
Hence, we assume 6(a) is a solution of a sixth degree equation, and look for its coefficients 
as polynomials in a. 



physical cut 




physical sheet 



Fig. 5: Analytic structure of b(a) in the 3-states dilute Potts model. Jumps of 6(a) 
across its cuts are shown. 



For the general dilute 3-states Potts model, the degree of these polynomials becomes 
quite high (up to 9), and we shall only write the algebraic equation satisfied by a and 6 in 
the non-dilute (a = 0) case: 



(3 6 6 + 6 7 (1 - (3 a) b 5 + (l3/3 7 2 a 2 - 6(-l + 4 7 2 ) a + 7 (2/3 + 9 7 //3)) 6 4 
+ ( - 12/3 7 3 a 3 + 4 7 (-6 + 7 7 2 ) a 2 - 4^(-6 + 2/3 2 7 + 3 7 2 ) a + ■ ■ ■ ) 6 3 

+ (4/3 7 4 a 4 - 6 7 2 (-5 + 7 2 ) a 3 + ^(9 - 54 7 2 + 10/3 2 7 3 - 15 7 4 )a 2 + • ■ • ) 6 2 

+ ( - 4 7 3 (3 + 7 2 ) a 4 - ^(18 7 - 24 7 3 + 4/3 2 7 4 - 18 7 5 ) a 3 + • • • ) 6 

+ (4 7 4 a 5 + 2- (17 - 18 7 2 ) a 4 + •••)= 

P 



(4.5) 
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Note that it was not known before that the standard resolvent of the Potts model uJs(b) 
(or a(b), cf Eq. fl2.5Q) satisfies a polynomial equation. This a fifth degree equation in a, so 
that a(b) has five sheets. The critical point is easily determined: 



7c 2 = 3 + V47 a\ = ^V 3 



which is compatible with what was found in [16]. The singularity of the resolvents is 



6-6,~(a-a,) 5/6 

( 4 -6) 

a-a*~(&-6*) 6 / 5 

which corresponds to a c = 4/5 theory. 

If we introduce the dilution, the equations become rather complicated, though it is still 
possible to work with exact analytic expressions. We find a critical line, with singularities 
of the type Q4.6|) , which ends with a tricritical point: 



a T = 2.44405 .. . (3 T = 2.9536 .. . 7 T = 6.09718. 
at which the singularity becomes 

b - b Teg ~ (a - a*) 7/6 

that is a c = 6/7 theory. 



(4.7) 



5. Q = 4 dilute Potts model 

For the Q = 4 case, it is easy to show that b(a) has an infinite number of sheets; of 
course, the method used above for Q < 4 then fails. In order to solve the model, the easiest 
procedure is to follow jl2| in which a similar "double saddle point equation" system was 
solved. One introduces an auxiliary function 

D(a) = 2b(a)-u A {a)-jj 

which, for Q = 4 (and only Q = 4), satisfies a two-cut Riemann-Hilbert problem: 

27 

D(a + iO) + D(a — iO) = —aa 2 + 70 — a E [ai, 02] lN 

p (5.1) 

D(a + iO) + D(a - iO) = a e [a 3 , +00] 
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where [ai, 02] is the physical cut and [03, +00] is the semi-infinite cut. D(a) can therefore 
be expressed via an elliptic parametrization in terms of O-functions. The critical line is 
expected to be found (just as in |12|) when the two cuts meet (02 = 03), that is in the 
trigonometric limit of the elliptic functions. The solution of ( |5.1| ) becomes then much 
simpler, and is of the form: 



D(a) = (a — 02) ( (a — ao) arctan * / — + a/ (a2 — ai)(ai — a) ) 

tv \ \ a± — a J 

where ai, a^-, clq and the relation defining the critical line are given by the asymptotics of 
D(a) as a — > 00 and the condition — a(a — 02) (a — ao) = — aaf> + 702 _ -g-- 

It is easy to see that positivity of the density of eigenvalues implies that ao > 02- In 
particular for a = (ao = +00), one finds 



7c 2 = 4(1 + v/TT^Ts) /3 C 2 = ^ 7c - 3 



More generally, along the critical line one has the strict inequality ao > a2 and the resolvent 
has the singularity 

b - 6re g ~ (a - a*) log 2 (a - a*) (5.2) 

or after inversion 

a - a reg ~ —2- — — (5.3) 
log (6-6*) 

On the other hand, at the tricritical point 

a T = 1.6523... (3 T = 2.42087... 7 T = 5.4602 .. . 
one has ao = 02, and this reflects in a change of the leading singularity: 

6 - 6 reg ~ (a - a*) log (a - a*) (5.4) 

or 

a — a reg ~ - — j-. (5.5) 

reg log(6-6*) v ; 



The behavior (|5.2|) and ( ^.4|) of the resolvent 6(a) is the usual one and is also found in other 



c = 1 models like the 0(2) model ]T5| , |2Ti[ ] ; however, the loop functions ( |5.3| ) (in particular 



the standard resolvent a(6) of the non-dilute 4-states Potts model) and ( p.5| ) are more 
unusual, and should be compared with the similar result found in [fT2"| . 
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6. Conclusion 

We have discussed a new method for solving multi- matrix models in the large iV 
limit, which relies on saddle point analysis and functional inversion relations satisfied by 
the unknown functions appearing in the saddle point equations. This method gives a 
new way of deriving equations for the resolvents of the model, and is a (much simpler) 
alternative to the method of loop equations for the determination of the disk amplitudes. 
It is not limited, like the orthogonal polynomials, to linear chains of matrices, as we have 
shown by solving explicitly the (Q + l)-matrix model representing the dilute Potts model 
on random surfaces. 

Furthermore, we have seen that the functional inversion can be interpreted here as a 
duality of the theory. In the Ising model this is simply the Z2 symmetry of the model; 
but more generally, in the dilute Potts model, the functional inversion relates different 
boundary conditions for surfaces with the topology of a disk. On the critical line, it 
exchanges two different loop functions with different critical exponents (g — > 1/g, while 
the bulk theory is unchanged, cf ( |1.7|) ) - similarly to the duality that exchanges Dirichlet 
and Neumann boundary conditions in open string theory. At the Q < 4 tricritical points, 
the universal part of the loop function turned out to be self-dual, whereas logarithmic 
corrections spoil the self-duality at Q = 4. 

One can hope that the ideas presented here will be applicable to various other problems 
of statistical mechanics on random surfaces, combinatorics, and asymptotics of orthogonal 
polynomials. 

Notes added 

After this work was completed, I became aware of the fact the critical point with 
singularities given by Eq. (|3.3|) is also the critical point of a percolation model, as noticed 
first in [|23| (I thank D. Jonhnston and M. Stathakopoulos for pointing this out to me). 
Also, it was claimed in a recent preprint []24| that using loop equations one could reproduce 
the polynomial equation ( fi~5| ) of the 3-states Potts model. 
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